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Model-independent determination of sharp features inside a star from 
its oscillation frequencies 
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It has been established earlier that sharp features like the base of the convective zone or the second helium ionisation 
zone inside a star give rise to sinusoidal oscillations in the frequencies of pulsation. The acoustic depth of such features 
can be estimated from this oscillatory signal in the frequencies. We apply this technique for the CoRoT"** frequencies of 
the solar-type star HD 49933. This is the first time that such analysis has been done of seismic data for any star other than 
the Sun. We are able to determine the acoustic depth of both the base of the convective zone and the Hell ionisation zone 
of HD 49933 within 10% error from the second differences of the frequencies. The locations of these layers using this 
technique is in agreement with the current seismic models of HD 49933. 
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1 Introduction 

A discontinuity in the derivatives of the sound speed in the 
stellar interior introduces a characteristic oscillatory signa- 
ture in the fr equencies of low degree modes of the form 
dGoughll 19901) 



sin(47rt/r + 0) 

where r is the acoustic depth of the layer, lying at a radius 
r, given by 

T — dr/c 

and 1/ is the frequency of oscillation. is a phase factor 

In solar-type stars, the base of the convective envelope 
and the Hell ionisation zones are the most prominent layers 
of such discontinuity. Each of these layers contributes to a 
sinusoidal variation in the frequencies of pulsation with the 
periodicity coiTesponding to the respective acoustic depth 
of the layer. The resulting combined signal in the frequen- 
cies can be fitted to a suitable functional form to extract the 
acoustic depths of these layers separately. Here we apply 
this technique to estimate the acoustic depths of these re- 
gions of the solar-type star HD 499 33 from its freque ncies 
obtained from the CoRoT satellite ( Baglin et al] |2006'). Al- 
though this m ethod has been used with great success f or the 



Sun (see e.g.. iBasu et alj|1994 iMonteiro et all 1 1994ft . this 



is the first instance of an application to a distant star. This 
was one of the main objectives of the CoRoT seismology 
programme dMichel et alj|2006l) . 



* anwesh@tifr.res. in 
** The CoRoT space mission was developed and is operated by the 
French space agency CNES, with the participation of ESA's RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. 



The significance of this method is that it does not de- 
pend on theoretical stellar models - the position of the sharp 
features are estimated directly from the frequencies them- 
selves. In fact, such independent determination of these lay- 
ers can help in the modelling of the star, as illustrated by 
Mazumdar (2005). We describe the technique used in this 
work and the results in the next two sections. 

2 Technique 

This oscillatory signal in the frequencies is quite small and 
is embedded in the frequencies together with a smooth vari- 
ation of the frequencies arising from the regular variation 
of the sound speed in the stellar interior. It can be enhanced 
by using the second differences (see, e.g., lBasu et al.lll994 , 
20041: iMazumdar & Antiall200lh , 



d^v{nj) = v{n - - 2v{n,t) + v{n + (1) 
instead of the frequencies themselves. Using the second dif- 
ferences serves the dual purpose of removing the smooth 
component in frequencies, and magnifying the oscillatory 
components. However, the errors in the frequencies are also 
enhanced in taking the differences, which might obscure the 
oscillatory signal. This is the reason why one cannot use 
even higher differences. 

The acoustic depths of the base of the convective zone 
(BCZ) and the second helium ionisation zone (HIZ), Tbcz 
and Thiz, respectively, can be obtained from the data by fit- 
ting the second differences to a suitable fun ction represent- 
ing th e oscillatory signals from these layers (iMazumdar & Antial 
200 lb. W e choose the following function along the lines of 
Basu et al. (2004): 



ao + 6o sin(47rz^TBcz + 0Bcz) 

+ {c(3+ci/v) sin(47ri^rHiz + 0Hiz) 



(2) 



®W]LEY 

JiTiterScieTice* 



© 0000 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



2 



Mazumdar, A. & Michel, E.: Acoustic depths of sharp features inside a star 



This function has three components corresponding to a 
residual smooth variation, and the two oscillatory compo- 
nents corresponding to BCZ and HIZ. W e do not use all the 
parameters in the function suggested by iBasu et al. (2004) 
to accurately represent the second differences because one 
cannot afford to use too many free parameters for a rel- 
atively small data set. The number of free parameters is 
optimised to strike a balance between a fair representation 
of the oscillatory signal and a reasonable x^- We find that 
the smooth component is fairly constant over the range of 
frequencies that we use, as is the amplitude of the signal 
from the BCZ. However, the amplitude of the signal from 
the HIZ varies more sharply with frequency, and thus re- 
quires at least on e frequency-dependent term. We note that 
Basu et alJ (120041) have shown that the exact form of the am- 
plitudes of the oscillatory signal does not affect the results 
significantly. 

The fit is carried out through a nonlinear ^ minimisa- 
tion, weighted by the errors in the data. The errors in the 
second differences are correlated, and this is taken care of 
by defining the using a covariance matrix. The contribu- 
tion of the errors are considered by producing 1000 realisa- 
tions of the data, where the mean values of the frequencies 
are skewed by random errors corresponding to a normal dis- 
tribution. The successful convergence of such a non-linear 
fitting procedure is somewhat dependent on the choice of 
reasonable initial guesses. To remove the effect of initial 
guesses affecting the final fitted parameters, we carry out 
the fit for multiple combinations of starting values. For each 
realisation, 100 combinations of initial guesses are tried for 
fitting the function above and the one with the minimum 
is accepted as the fit for that particular realisation. 

The median value of each parameter for its distribution 
over 1000 realisations is taken as the average. The error in 
the parameter is estimated from the range of values cover- 
ing 68% area about the median (corresponding to Icr error). 
As an example, the histogram of distribution of the param- 
eters Tbcz and Thiz over 1000 realisations are illustrated in 
Figs. |2] and [3] The quoted errors in these parameters reflect 
the width of these histograms on two sides of the median 
value. 

3 Results 
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Frequency v (mHz) 

Fig. 1 Second differences of frequencies of HP 49 933 
from 180 days of CoRoT data ( Benomar et al. 20091) are 



plotted against frequencies for ^ = {blue points), £ = 1 
{magenta points), and a few t — 2 {red points) modes. The 
fit of the oscillatory function is shown by the brown solid 
line. The cyan dotted line shows the fitted oscillatory signal 
due to the base of the convection zone (BCZ) and the green 
dashed line shows the signal due to the second helium ioni- 
sation zone (HIZ) 
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Fig. 2 Histograms of fitted values of acoustic depths, Tbcz 
{red) and Thiz {blue), obtained for 1000 different rea lisa- 
tions of the data according to iBenomar et al. ( 20091) are 
shown. The median values are shown with the dotted lines. 



Our primary data set is that of iBenomar et al I (l2009l). To 

check for consistency we also use the data set of lOruberbauer et\ail[i es set equal to the median values obtained from 1000 
(|2009h. and the th eoretical frequencies o f a model of HP 49933 realisations. The histograms of Tbcz and Thiz values over 



Peheuvelsll2010l: see also lGoupilll2010l) . 

The fit of the fun ction given in Eq.lHis shown in Fig.[T] 



This set con- 



for the data set from IBenomar et al. 
sists of 31 second differences computed from the individual 
frequencies. We did not use frequencies which had errors 
of more than 3 /^Hz. The data represented in this figure cor- 
responds to the mean frequencies. In reality, different re- 
alisations of the data are fitted independently. The figure 
also shows the function given in Eq. |2] with the parameter 



all realisations are shown in Fig.|2] The median values are 
indicated, and the errors are estimated by considering the 
extent of the histogram around the median such that 34% 
of the cases are covered on each side. Since the histograms 
are asymmetric, naturally this yields unequal error bars. We 
consider these to be representative of Icr error on the me- 
dian value of the parameter concerned. The average am- 
plitude of the oscillatory components corresponding to the 
two regions, A^cz and A^iz, are evaluated in terms of the 
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Fig. 3 Similar plots as Figs. [T]and|2l but with the frequencies taken from the 60 day CoRoT data (IGruberbauer et al 
20C)9l) . 



freq uencies themselves, inste ad of the second differences 
(see Mazumdar & Antial2001 , for details). These results are 
summarised in Table [U 

To check for consisten cy, we also used the data set from 
IGruberbauer et al. ( 20091) . which yielded 20 second differ- 
ences. This data set yields very similar results, shown in 
Fig.[3]and listed in Table[T] although with higher error bars, 
showing that the method works even with lesser data. 

Lastly, we also fit Eq. |2] to the theoretical frequencies 
computed from a s tellar model of HD 49933 (IDeheuvels 
20101: iGoupifcoiOh . The actual model values of the acous- 
tic depths, Tbcz and Thiz, are remarkably close to the val- 
ues obtained from this fit (see Table [T]|. This is not entirely 
unexpected, since the stellar model is obtained by match- 
ing the seismic frequencies themselves. However, the meth- 
ods used in the seismic modelling do not explicitly involve 
determining the depths of the base of the convective enve- 
lope or the second helium ionisation region. In contrast, our 
method focuses directly on these two layers, irrespective of 
the str atification in the re st of the star. However, as pointed 
out by 'Mazumdai ( 2005 ). the location of these layers can- 



not be independent of the structure of the rest of the star, 
and is thus linked implicitly to the detailed model. But the 
process of fitting the acoustic depths does not involve stellar 
modelling and thus the results obtained in this method are 
completely model-independent. 

4 Summary 

We have determined the acoustic positions of the base of the 
convective envelope and the second helium ionisation zone 
in the star HD 49933 from the oscillatory signal in its fre- 
quencies observed with CoRoT. The method involves fitting 
a suitable function to represent the two oscillatory signals 
arising from the two zones to the second differences of the 
frequencies. We find almost similar results for two data sets, 
one corresponding to only the initial 60-day run of CoRoT 



a. 




Frequency f (mHz) 



Fig. 4 Second differ ences of frequencies of a theore tical 
model of HD 49933 jDeheuvelsl l201ot lGoupi]|l201(jl) are 
plotted against frequencies for i = 0, £ = 1, and £ = 2 
modes. The composite fit of the the oscillatory function, as 
well as its two components are shown, similar to Fig. [l] 



and another a 180-day combined time-series obtained with 
a long run of CoRoT. This is the first instance where such a 
technique has been successfully applied to a star other than 
the Sun. 

Although these results are independent of any stellar 
modelling, they match well with a seismic model of the star 
In fact the acoustic depths of such layers of sharp changes 
in the sound speed inside t he star extracted in this method 
can help to model the star (lMazumdaj2005h . 

The same technique can be used to estimate the acoustic 
depths of such layers in other stars. We have established that 
a time series of moderate length, as obtained by CoRoT, is 
sufficient for the successful application of this technique. A 
longer time series, as expected from the Kepler mission, for 
example, would further reduce the errors on the extracted 
acoustic depths. 
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Table 1 Results of the fitting of second differences to Eq.|2]for two data sets and a theoretical model. The quantities in 
parentheses represent the errors in each parameter for the first two rows, and the actual theoretical model values in the third 
row. 



Data Set 


No. of 


Acoustic Depth 


Average Amplitude in v 






data pts. 


Tbcz (s) Thiz (s) 


^bcz(mHz) Ahiz(mHz) 




CoRoT 180 d fBenomar et al. 2009) 


31 


2150 (+182/-257) 836 (+56/-52) 


0.57 (±0.30) 1.29 (±0.30) 


1.03 


CoRoT 60 d (Gruberbauer et al. 2009") 


20 


2253 (+2897-226) 845 (+191/- 142) 


0.52 (±0.20) 0.93 (±0.32) 


0.83 


Model (Deheuvels 2010: Goupil 2010) 


57 


1995 (1999) 826 (883) 


0.46 1.35 
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